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1 INTRODUCTION 



Since the first dete ctions in thermal infra red dAumann et al ] fl984h 
and scattered light dSmith & Terrilel 19841) . it has become clear that 
many main-sequence stars are surrounded by discs of dust grains. 
The dust grains have only a short lifetime compared to the age of 
the star, being ground down in collisional processes until they are 
small enough to be ejected from the system by radiation pressure 
jDominik & Decinl2003l ; IWvattl2005h . The existence of a dusty de- 
bris disc therefore implies a large reservior of planetesimals, par- 
ent bodies with longer collisional lifetimes whose collisions are 
producing the observed dust. Such planetesimals must be collid- 
ing with sufficient relative velocity to produce enough dust to be 
observable, and so the disc must have some degree of dynami- 
cal excitation, with non-zero eccentricities and inclinations caus- 
ing high relative velocities. Although the formation of planetesi- 
mals is still not understood, it is usually assumed that they form 
on coplanar, near-circular orbits in the protoplanetary disc (e.g., 
Kenvon & Bromlev 2008). Therefore, there must be some means 
of exciting the planetesimals' eccentricities for the disc to become 
visible. This is referred to as stirring. 

The origin of debris disc stirring is not yet known and has 
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ABSTRACT 

Detectable debris discs are thought to require dynamical excitation ('stirring'), so that plan- 
etesimal collisions release large quantities of dust. We investigate the effects of the secular 
perturbations of a planet, which may lie at a significant distance from the planetesimal disc, to 
see if these perturbations can stir the disc, and if so over what time-scale. The secular perturba- 
tions cause orbits at different semi-major axes to precess at different rates, and after some time 

Across initially non-intersecting orbits begin to cross. We show that t CIOSS oc <^J^ C / {rn^e^a^, 
where m v \, e p \ and a p \ are the mass, eccentricity, and semi-major axis of the planet, and adisc 
is the semi-major axis of the disc. This time-scale can be faster than that for the growth of 
planetesimals to Pluto's size within the outer disc. We also calculate the magnitude of the 
relative velocities induced amongst planetesimals and infer that a planet's perturbations can 
typically cause destructive collisions out to 100's of AU. Recently formed planets can thus 
have a significant impact on planet formation in the outer disc which may be curtailed by 
the formation of giant planets much closer to the star. The presence of an observed debris 
disc does not require the presence of Pluto-sized objects within it, since it can also have been 
stirred by a planet not in the disc. For the star e Eridani, we find that the known radial velocity 
planet can excite the planetesimal belt at 60 AU sufficiently to cause destructive collisions of 
bodies up to 100 km in size, on a time-scale of 40 Myr. 

Key words: circumstellar matter - planetary systems: formation - planetary systems: proto- 
planetary discs - stars: individual: e Eridani - stars: individual: Fomalhaut. 



yet to be thoroughly investig ated. The most comprehen sive model 
of this is self — stirring dKenvon & Bromlev f l2008t and refer- 
ences therein) where a planetesimal disc evolves due to mutual low- 
velocity collisions resulting in the growth of the planetesimals. In 
this model the disc is stirred, and collisions become destructive, 
when the largest planetesimals reach Pluto's size, and their gravita- 
tional perturbations excite the remaining smaller bodies. 



Such planet formation models, however, ignore the effect of 
the formation of massive planets on the evolution of material. Yet 
we know that large numbers of massive planets exist, including 
withi n debris disc systems (e.g., Benedict et al. 2006: iKalas et al.l 
2008). The existence of a massive planet in a disc system can 
have several consequ ences, including migration of the planet 
dFernandez&ldll984h and scattering of embryos from the inner 
disc ne ar the planet to the outer disc (Goldreich. Lithwick & Saril 
l2004al) . The formation of a planet may also be able to speed planet 
formation exteri or to the planet when the gas disc is still present 
dThommesl2 005). In this paper, we consider the effects of a planet's 
gravitational secular perturbations on the disc, showing that they 
can stir the disc on timescales of typically several 10s of Myr. 

This paper is organised as follows. In §2 we review the secular 
theory we are using to model the effect of the planet on the disc. 
In §3 we derive an analytical estimate for the time the planet takes 
to stir the disc. In §4 we investigate the relative velocity distribu- 
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tion imposed by the planet. In §5 we discuss the model we use to 
ascertain the outcomes of collisions. In §6 we discuss the implica- 
tions of planet stirring for debris disc evolution. Finally, in §7 we 
summarise the paper. 



2 DYNAMICS 

We consider the orbital evolution of planetesimals, treated as mass- 
less test particles, under the perturbing influence of N massive 
planets. We shall consider only secular perturbations. In the fol- 
lowing, unsubscripted orbital elements refer to test particles and 
subscripted to the planet(s). 

We use classical Laplace-Lagrange secular theory to model 
the long-term orbital evolution. This is valid if planetesimals are 
not near mean motion resonance and eccentricities are low. The 
planetesimal's complex eccentricity z = e exp (izu) can be decom- 
posed into a forced component zi and a proper component z p : 

z(t) = z f (t)+z p (t) (1) 



E 



.4 



9i 



■ exp [i (gd + /3i)] + e p exp [i (At + P)] 



dWvatt et alj 19991) . Here, and eji are the eigenvalues and eigen- 
vector components of the solution for the planets, independent of 
the planetesimal's location, and Pi, P, and e p are constants of in- 
tegration; e p is the proper eccentricity. The constants Aj are given 
by 
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(2) 



with ?7i j the mass of the j-th planet, and A is given by A 
Ef =1 B„with 

B ^^va 



+VW^^a 3 a i6 «(a 3 .)/4, 



(3) 



ay = a.j/a for an interior planet and ay = a/aj for an exterior 
plane t, (a) are the Laplace coefficients (see lMurrav & Dermottl 
1999), and ay = 1 for an interior planet and ay = ay for an 
exterior planet. 

The behaviour of the complex eccentricity is to precess at a 
rate A on a circle of radius e p about the forced complex eccentricity 
Zf. For more than one planet, Zf is itself evolving in time. If there 
is only one planet, its orbit is fixed, and so is Z{. 

Henceforth we shall concentrate on the case of a single planet. 
We shall return briefly to the multi-planet case in £|6.4.3| Denoting 
the planet's elements now by p i, we have 



z(t) 



"3/2 

b^ 

°3/2 



e p i exp (ixvpi) + e p exp [i (At + 0)] . 



(4) 



We shall later concentrate on the case where the planetesimals' or- 
bits are initially circular. In this case, e p = ef and P — —za p \. 

U sing the leading-order app roximations to the Laplace coeffi- 
cients ([ Murray & Dermottlll999l) : 
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Figure 1. Secular precession time-scale t scc = 2tt/A (top) and forced 
eccentricity ef (bottom) for a planetesimal on an initially circular orbit at 
a = 15 AU, with a p j = 5AU, m* = IMq, and ji = m D i/m <. and e D i 
varyin g as shown. Symbols: N-body simulations using RADAU (Everhart 
1985). Solid line: analytical prediction from Laplace-Lagrange theory. Dot- 
ted line: analytical prediction from Laplace-Lagrange theory with Laplace 
coefficients approximated to leading order in a (Equations f5J. Dashed line: 
analytical prediction from Heppenheimer-Kaula theory. 



in the limit of small a. 

In Figure Q} we compare t BCC as predicted by Laplace- 
Lagrange theorytothat obtained by numerical integration (using 
RADAU, Everhar jl 1985b . Precession timescale t scc was calculated 
by fitting a sinusoid to the eccentricity output from the integratiorQ. 



1 Orbital elements from the simulations are derived from canonical Jacobi 
coordinates (motion of planet referred to the star, motion of test particle 
referred to the star-planet barycentre). This eliminates variati ons in eccen - 
tricity (of order ifj,/^/a) on the Keplerian timescale I Lee & Pealel l2003l) . 
which exceed the forced secular eccentricity at large distances from the star. 
Use of Jacobi coordinates also ensures no dependence of secular behaviour 
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While Laplace-Lagrange theory is very accurate at small planetary 
eccentricities, it overestimates the precession time-scale and forced 
eccentricity at high e P i. For high eccentric ity planets, we therefore 
adopt the theory of Heppenheimen a 19781) based on the disturbing 
function expansion of Kaulal d 19621) . This is an expansion in a with 
no restriction on e v \, in contrast to the Laplace-Lagrange expan- 
sion which is for small e p i without restriction on a (provided there 
is no mean motion resonance). This predicts qualitatively the same 
behaviour, with forced eccentricity 



ef 



5<ae p i 



4 1 



-pi 



and precession rate 



A = y Qm i ,/a z - 
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From Figure (TJ, we see that, while the Heppenheimer-Kaula 
theory provides a better description of secular behaviour of t sec 
than Laplace-Lagrange theory for high e p i, it actually gives a worse 
fit for ef. We therefore omit the factor (l — e pl ) 1 from ef from 
now on, as an empirical correction. We also see in Figure Q} that 
the Heppenheimer-Kaula solution is overestimating t scc for low 
e p i. This is due to higher order terms in a which the Heppenheimer- 
Kaula solution does not take into account. For comparison, the 
Laplace-Lagrange solution is also plotted to leading order in a. 
As a becomes small, the Heppenheimer-Kaula solution performs 
better at small e p i. Note that at high e p i there is also a dependence 
of tsec and ef on /j, beyond that predicted by either theory; this is 
probably due to both theories being based on treating the perturba- 
tions only to leading order in fi. 

Figure ([2} shows the dependence of i scc on a for three dif- 
ferent planet masses. The Laplace-Lagrange theory describes well 
the behaviour of test particles, except for those near mean motion 
resonance with the planet. The effect of resonances covers a larger 
range of semi-major axes as planetary mass is increased. The anal- 
ysis we present below is based on a purely secular theory and so 
does not apply for planetesimals near resonance. For /i = 10 -3 
and a p i = 5 AU, this excludes planetesimals within 8.5 AU and 
those and those at 10-11 AU. 

In subsequent sections we shall derive expressions for the rela- 
tive velocities in the disc and the time-scale for disc stirring correct 
to first order in a. Henceforth we use Equation ^ for the preces- 
sion rate, and Equation Q for the forced eccentricity. 

The secular precession rate A is a strong function of a (Equa- 
tion^ Figure O. Furthermore, the forced eccentricity ef also de- 
pends on a (Equation [6]l. The eccentricity evolution and apsidal 
precession of planetesimals at different semi-major axes is differ- 
ent, and this differential precession raises the possibility that the 
initially non-intersecting orbits may cross after a certain time t CIOSS . 
We investigate this possibility in the next section. 



3 TIME-SCALE FOR ORBIT CROSSING 

We are interested in whether a planet's secular perturbations can 
cause neighbouring planetesimals to collide. Consider two plan- 
etesimals at semi-major axes ai and a2, with a± < 02- The 



on initial mean longitude at the same semi-major axis, which is not the case 
for heliocentric coordinates. 
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Figure 2. Secular precession time-scale t SEC as a function of planetesi- 
mal's semi-major axis, with a p i = 5AU, m* = IMq, fi = 0.001, and 
e p i = 0.1. Solid line: L-L prediction. Crosses: simulation results. Note the 
presence of strong mean motion resonances at the indicated locations. Plan- 
etesimals which were ejected within the integration time (10 5 , 10 6 , and 
10 7 years for fj, = 10 — 2 , fj, = 10 -3 , and fi = 10 -4 respectively), mostly 
close to the 2:1 resonance with fi = 10 -2 , are not shown. The analysis in 
this paper is not valid for planetesimals in resonance. 



orbits are interecting at a time t if n (ai, ei(t),U7i(t),6) ^ 
Ti (02, e.2{i), W2(t), 9) for some true longitude 8. Writing the 
equation of an ellipse to first order in eccentricity, this condition 
becomes 



ai [1 — ei cos(cci — 9)] ^ a,2 [1 — e2 cos(ci72 — 9)] 



(10) 



Now, we consider closely separated orbits, so that 02 = a\ + da. 
We then require 
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+ cii cos (uji - 
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de p AA 

— aie p it sin ( vox — o) — — 
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aie p ,i sin (zui 



(11) 



All the terms are 0(e) except for the first term and the term 
— aie p ,itsin (01 — 9) 4^ which becomes order unity for suffi- 
ciently large time t, provided that 4^ 7^ (this condition is sat- 
isfied for single-planet systems, but multi-planet systems can have 
4^- = 0; see 36.4.3t . If 4^ = then second-order terms in the 

da " ' aa 

Taylor series must be considered; the term of order unity will be 
oc e D idat 2 ^4-. 

^' da- 

We now use the leading-order expansions for the Laplace 
coefficients (Equation [5j, so with small eccentricity e and 
planet/planetesimal semi-major axis ratio a we have 



> 1 + bAte p sin (At + (3 - 9) + O (e) 



(12) 



for orbit intersection, where b — 7/2 for an internal perturber and 
b — —3/2 for an external perturber. The second term is of order 
unity when t is of order e p . 

At the crossing time t — t CTOSS , 



1 + bAt cross e p sin (At cross + (3 - 9) + O (e p 



0. 



(13) 



4 A. J. Mustill & M. C. Wyatt 



NOW, | Sin (Across + P ■ 

for some 6, so 
1 



\b\Ae 



and it will attain its maximum 



(14) 



Note that this has no explicit dependence on ef, but can have a 
dependence on ef if e p depends on ef. In the remainder of this 
paper, unless otherwise stated, we assume the orbits are initially 
circular; this means that ef = e p . Then, for an internal perturber, 
we have 
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while for an external perturber, we have 

\3/2 
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(16) 



Orbit crossing begins more quickly if the perturbing planet is 
more massive or eccentric, and the time-scale is a strong function 
of planetesimal semi-major axis. 

Equation 116t i s in good agreement with equation (14) of 
Theb ault et ail (12006 ). which was derived as an empirical law based 
on a slightly modified secular theory and N-body simulations, for 
test particles on circumprimary orbits in a binary system (i.e., ex- 
ternal very massive perturber). The chief difference between their 
result and ours is for high eccentricity. 



4 RELATIVE VELOCITIES 

Now that we have established the time-scale on which orbits cross, 
we proceed to examine the relative velocities of planetesimals un- 
dergoing collisions. We shall examine how the distribution of rela- 
tive velocities evolves with time. 

For planetesimals with randomised, uniformly distributed ap- 
sides, we might expect the mean relative velocity in a collision to 
be given by 

(v Icl ) = c (e) u kcp , (17) 

where c is a constant of order unitjQ. The value of c depends on 
the specific definition of relative velo city being used, and the un- 
derlying eccentricity distribution (see lLissauer & Stewart! 1993h ; in 
our case we wish to know the mean velocity of a planetesimal 
relative to others in the swarm, for which c = \/5/4 when the 
planetesimals' eccentricity follows a Rayleigh distribution, which 
arise s from the mutual g ravitational scattering of planetesimals 
(e.g., flda & Ma Mnolll992l) . However, the eccentricity distribution 
arising from a planet's secular perturbations cannot be assumed to 
be Rayleigh, because the physical process exciting the eccentrici- 
ties is very different (long-range secular perturbations vs. mutual 
gravitational scattering). Furthermore, the apsides are constrained 
by zu 6 (zu p \ — t;-k, vo-pi + ^tv). This is because, with the or- 
bits initially circular, the complex eccentricity starts at the origin 
of the complex plane and precesses in a circle around z\. Thus it 



For inclined orbits the inclination also plays a role. 




Figure 3. Maximum (solid line), mean (dotted) and median (dashed) rela- 
tive velocity experienced by a planetesimal, as a function of time. The plan- 
etesimal is at 15 AU, and the planet at 5 AU with eccentricity e p i = 0.1 and 
mass m p i = 10 — 3 Mq, orbiting a solar mass star. 



is restricted to the half-plane containing zf. Finally, Equation d 1 7b 
only applies locally, whereas in reality a planetesimal on an eccen- 
tric orbit can collide with others over a range of semi-major axes 
(see iThebault & Doressound iram 2003, for a discussion). In this 
section, we therefore examine the velocity distribution imposed by 
planetary secular perturbations, beginning with an estimate for the 
range of semi-major axes over which collisions can occur. 

First, we estimate the maximum radial excursions of planetesi- 
mals evolving under secular perturbations. Consider Planetesimal 1 
located at semi-major axis a%, with eccentricity ei = |aie p i/a p i, 
the maximum attainable under the planet's secular perturbations 
starting from an initially circular orbit. We wish to find the greatest 
semi-major axis of an exterior planetesimal, Planetesimal 2, such 
that Planetesimal 2's orbit can intersect that of Planetesimal 1. If 
longitudes of periapse could take any angle, this would occur when 
Planetesimal 1 was at apapse and Planetesimal 2 at periapse, with 
the orbits tangent, and the apsides antialigned.However, the sec- 
ular solution also imposes a restriction on the longitude of peri- 
apse zu, restricting it to the range zu G (zu p \ — tt/2, zu p i + 7r/2). 
Because of this, the maximum semi-major axis for Planetesimal 2 
must come when Planetesimal 2's orbit is at its lowest eccentricity, 
i.e., e2 = 0. Denoting the difference between the semi-major axes 
of the orbits by Aa, we find, to lowest order in e p i and a, that the 
maximum separation of intersecting orbits is 



(A«; 



max = —aiae p i, 



(18) 



which is simply the maximum radial excursion of Planetesimal 1. 

We now numerically calculate the distribution of relative ve- 
locities which a planetesimal at ai experiences, as a function of 
time, assuming that all the planetesimals evolve deterministically 
under the Laplace-Lagrange secular solution described in Sj2] start- 
ing on initially circular orbits. For interactions with planetesimals at 
different se mi-major axes, the rel ative velocity is calculated using 
formulae in lWhitmire e t al. ( 1998). Note that two confocal ellipses 
can have two intersection points; we calculate the relative velocity 
at each point and use both in the analysis below. 

Figure {3} shows the relative velocities encountered for a fidu- 
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cial case of a p i = 5AU, a± = 15 AU, e v \ = 0.1, = IMq, 
and m p i = 1O _3 M0, for i = to t = 10t cr oss- The surface 
density of planetesimals is assumed to be constant. We also calcu- 
lated the averages with a E oc a~ 3 ^ 2 profile with negligible dif- 
ference in the derived relative velocities. We see that planetesimals 
begin on non-intersecting orbits, and evolve for t ~ tcross before 
orbits begin intersecting. For these parameters, orbit-crossing be- 
gins slightly sooner than t = t C ross because tcross was derived in 
the limit a<l, and here we have a = 1/3. We also see that the 
maximum relative velocity experienced by a planetesimal rises to 
« 2ef«kc P . This is the relative velocity of a planetesimal with ec - 
centricity 2ef relative to a circular orbit l lLissauer & Stewart! 1993b . 
but can be achieved for other configurations too. The average rel- 
ative velocity - mean or median - is m 1.4ef«ke P - Therefore, we 
can in fact use Equation dl7| (, with c = 2 if we wish to use the 
maximum relative velocity, and c = 1.4 if we wish to use the aver- 
age. We also note that, while the analytical orbit-crossing criterion 
only guarantees intersection of infinitesimally separated orbits, we 
found that soon after orbit crossing began, the planetesimal's orbit 
intersected the orbits of planetesimals over the whole range permit- 
ted by Equation d 1 8 b - 

We can now derive expressions for the maximum relative ve- 
locity imposed on a disc by a planet. Combining Equations l[6} 
and dl7t with the expression for the Keplerian velocity v^ cp — 
\J Qrrit/a, we get the following expressions for the maximum rel- 
ative velocities: 
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for an internal perturber, and 
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(20) 



for an external perturber. If it is desired to work with the mean 
or median relative velocity, these numbers should be reduced by a 
factor of m 0.7. 

In Figure {4]l we compare the relative velocity distribution im- 
posed by a planet's secular perturbations with that arising from a 
Rayleigh distribution of eccentricities. We see a clear difference 
in that the Rayleigh distribution gives a tail of high- velocity col- 
lisions, in contrast to the planet-stirred distribution's well-defined 
maximum. In the remainder of this paper, we shall use the maxi- 
mum relative velocity. 



5 COLLISION OUTCOMES 

We now need to relate the relative velocities in the disc to the out- 
comes of collisions between planetesimals. When relative veloci- 
ties are low, collisions between planetesimals result in net growth 
of the larger body. When relative velocities are high, collisions re- 
sult in the bodies fragmenting. In the former regime, smaller bod- 
ies build up to eventually produce planetary embryos. In the lat- 
ter regime, planetesimals are ground down in a collisional cascade 
with the production of large quantities of dust which can be ob- 
served as a debris disc. We now proceed to quantify the minimum 
velocities needed for these erosive collisions to occur, and to find 




v rel/ e f v kap 



Figure 4. Black: cumulative frequency distribution for relative velocities in 
a planet-stirred population at t = lt G ross (solid line), t = 2t cross (dot- 
ted), and t = 10t crO ss (dashed). Grey: cumulative frequency distribution 
for relative velocities in a population with a Rayleigh distribution of eccen- 
tricities. The mean eccentricity for the Rayleigh distribution is the same as 
that for the planet-stirred case at t = 10t C ross- 



which combinations of planetary parameters can lead to such ve- 
locities. 

Consider two bodies colliding. The outcome of the collision 
depends on a large number of factors including the bod ies' compo- 
sition, shape etc. (see, e.g jrlousen & Holsapplefll990l) . but for our 
purposes the most important considerations are the relative velocity 
and the kinetic energy of the impact. 

We are interested in the transition from net accretion to net 
erosion, and so use the condition Qr > Qrd, where 



0.095 J ke 



+0.025 J kg' 
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(21) 



for weak aggregates dStewart & Leinhardtll2009l) . Here, Rci is the 
radius of a sphere with mass equal to the total of the two bod- 
ies' masses mi + m2, and a density of lgcm" 3 , and Qr, = 
0.5ur Gl mim2/(rai +m2) 2 is the reduced mass specific kinetic en- 
ergy. We use the results for weak aggregates rather than strong rock 
because firstly we are primarily interested in outer discs which have 
a high proportion of ices, and secondly we are considering bodies 
formed by sequential accretion of small er ones, which results in 
porous aggregates rather than monoliths dWeidenschillindl 19971) . 

Here we define the catastrophic disruption threshold velocity, 
v* e \, to be the velocity above which collisions between equal-sized 
bodies no longer result in one of these bodies gaining mass. It is 
given by 



0.8 



R 
80 m 



+ 0.2 



R 



-\(22) 



where we have converted Rci into a physical radius R for equal- 
sized bodies at a density of 1 g cm -3 . This function has a minimum 
of w 1ms" 1 atRm 80m. 
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Equation J22b gives a lower limit to the relative velocity 
needed to destroy a body of radius R. Because the velocity dis- 
tribution excited by a planet's secular perturbations has a definite 
maximum, which decreases with an outer planetesimal belt's semi- 
major axis, we see that a planet's secular perturbations will be un- 
able to cause catastrophic collisions in a disc beyond some critical 
semi-major axis 
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for an internal perturber, and within a critical semi-major axis 
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( fl pi 



V10AU 



VIM© 
<d(-R) 

lkms" 1 



(24) 



for an external perturber. 

We can also estimate the size of the largest body which a given 
planet can destroy at a given disc radius a. For R 3> 80 m, we 
approximate Equation J22b by w* e i ~ 0.26(i?/80 m) ms _1 ; com- 
bining this with Equation d 1 9b for the maximum relative velocity 
gives 
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74km (2SL 
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X (lOAu) 
for an internal perturber. 
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6 DISCUSSION 

6.1 Ability of a planet to stir a disc 

We can now discuss the effects of planetary secular perturbations 
on an initially dynamically cold planetesimal disc. We have seen 
how the perturbations induce orbit-crossing and excite relative 
velocities between planetesimals (SQ. Whether this is sufficient to 
stir a disc depends on the typical size of bodies in the disc (see 
Equations 1231 and |24l . If the planetesimals in the disc have grown 
to ~ 80 m in size, then a planetary perturber can stir the disc at 
radii 

1/3 

a < a* = 1800 AU eT [ ) ( ) . (26) 



Because the 80 m bodies are the weakest, Equation d26| l determines 
the greatest range of a planet's secular perturbations for disc stir- 
ring. If R S> 80 m, we can write 



a* = 170 AU e 2 p { 3 



R 



10 km 



-2/3 



i in 



1/3 



Opl 

1 AU 



2/3 



(27) 



We illustrate Equation d26t in Figure QJ. This figure shows 
the maximum semi-major axis a* at which 80 m bodies can be de- 
stroyed as a function of a planet's position in a-e space. We see that 
all but very low eccentricity planets close to their star have the po- 
tential to stir discs out to at least 30 AU. Hot Jupiters are therefore 
ruled out as potential disc stirrers. 

We note that this constraint is independent of the planet mass. 
The planet mass affects the time for perturbations to act, but not 
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Figure 5. Maximum semi-major axis a* at which 80 m bodies can be de- 
stroyed, as a function of planetary semi-major axis and eccentricity (Equa- 
tion ^) . This is independent of planet mass. Stellar mass = 2Mq. 



their eventual effects, assuming that the planet's secular perturba- 
tions are the only source of dynamical evolution. For low mass 
planets, the secular timescale may be comparable to the timescale 
for col lisional damping of eccentricity ( Goldreic h. Lithwick & Saril 
l2004bl) . 

6.2 Comparison to self-stirring models 

In the self-stirred model (see §1), the disc begins in an unexcited 
(e ~ 10~ 4 ) state, composed of sub-km planetesimals. These grow 
through collisions until they reach sizes similar to Pluto's, at which 
point their gravitational perturbations stir the disc, causing the disc 
to brighten. On the assumption that both the self-stirring and planet- 
stirring models are accurately describing the behaviour of the disc, 
a key question for disc evolution is which occurs sooner. In this 
section we compare the time-scales for the two processes to occur. 

The time taken to form Pluto-sized bodies at a given radius 
through core accretion, tp\, is proportional to the orbital time-scale 
divided by the disc surface density at that radius. Base d on exten- 
sive numerical simulations, Ken von & Bromlevl J2008h find 

3/2 



tpi = 145a; n 



( a \ 3 /2M y 
V 80 AU / V m * J 



Myr. 



(28) 



Here, x m parametrises the disc surface density in such a way as 
to account for the propensity of more massive stars to have more 
massive discs: the surface density of disc solids is given by 

E = Eox m (m,/M Q )(a/a )~ 3/2 , (29) 

where Eo = 0.18gcm -2 corresponds roughly to the minimum 
mass solar nebula density at ao = 30 AU. We assume that if gi- 
ant planets are present then the planetesimal surface density is not 
depleted below the primordial value by processes such as planet- 
planet scattering or planetesimal-driven migration. While depletion 
would lead to longer self-stirring timescales, such processes may 
also excite the eccentricities of remaining planetesimals, efficiently 
stirring the disc at earlier times. Here we concentrate exclusively 
on the planet's secular perturbations. 
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Figure 6. Comparison of the time-scales for stirring in a system with a 
planet at 5 AU with eccentricity 0.1, a 2 solar-mass star, and an extended 
planetesimal disc beyond 10 AU. Solid lines show planet-stirring time- 
scales for rvipi = 10~ 4 , 10 — 3 , 10 -2 Mq. Dotted lines show self-stirring 
time-scales tpi for disc masses given by x m = 0.1, 1, 10. More massive 
planets stir the disc more quickly; more massive discs stir themselves more 
quickly. 



Figure 7. Illustration of constraint 1301 . that if planet-stirring is competing 
with self-stirring then the disc is planet-stirred precisely when the time- 
scale for planet stirring is shorter than that for self stirring. Lines are plotted 



at fixed a 



pi 



5 AU, 



3 p1 



0.1, 



2 Mq , and m^i as shown. At a 



'pi 



given disc radius, the planet can stir the disc for disc surface densities on 
the lower (shaded) side of the line. 



Figure l[6) compares the self-stirring time-scale ipj given by 
Equation l |28t with the time-scale for planet stirring t CIOSS given 
by Equation d 1 5 b - At a given planetesimal belt semi-major axis, tpi 
depends only upon the disc surface density, while t C ro SS depends on 
planet mass, semi-major axis and eccentricity. Furthermore, both 
time-scales have a different dependence on stellar mass. As t C ross 
has a stronger dependence on a, for any given planet parameters 
and disc mass, there exists a radius beyond which the disc will stir 
itself before the planet can stir it. For example, for a Jupiter-like 
planet and a disc with x m = 1, t CIOSS < tpi out to about 13 AU, 
whereas for a heavier 10 Jupiter mass planet t cr oss < tpi out to 
around 60 AU. 

With this in mind, we define a new parameter $ which de- 
scribes the boundary between planet-stirred and self-stirred regions 
of the disc, assuming that planet-stirring and self-stirring are the 
only two mechanisms operating. The disc is planet stirred if 

a<$ = 630AU (l-e^-e^^) 273 

/ dpi \ 2 /m*\" 4/3 _ .77 

and self-stirred otherwise. 

This is illustrated in Figure Q, showing the disc surface den- 
sity below which a planet at 5 AU, with e p \ — 0.1, can stir faster 
than the disc can form Plutos. If the planet is of Jupiter mass, it can 
stir discs with Xm < 1 at 20 AU. 

If the disc is planet-stirred before it self-stirs, then there will 
be implications for planet formation. If the largest bodies that the 
planet can destroy (Equation [25} are larger than the largest bodies 
in the disc, then further growth of planetesimals will be difficult 
if not impossible. Even if the largest planetesimals are unable to 
be destroyed, their growth rates will be reduced as the increased 



velocity dispersion amongst the planetesimals reduces gravitational 
focusing factors, leading to longer collision time-scales. 

6.2.1 Observables of stirring models 

If the disc has an inner hole then we might expect a low level of 
dust production, hence IR luminosity, until the planet's secular per- 
turbations cause orbits to cross at the inner disc edge. The disc 
would then brighten (at a given wavelength), before getting dim- 
mer again at that wavelength as the region of peak dust produc- 
tion moves outwards. This could explain the observed incidence 
of excess IR emission and the fractional luminosity for young A- 
type stars, which both appa rently increase wi th age, peaking at 
10-20 Myr before declining dCurrie et alj[2 008). This behaviour is 
qualitatively similar to that predicted by self-stirring, where the re- 
gion of peak dust production moves outwards as Pluto-sized bodies 
form at progressively larger radii. The peak at 20 Myr cannot there- 
fore be taken as evidence for one particular type of delayed stirring 
over another. 

If most planetesimal discs are extended rather than being nar- 
row rings, we would expect to see an increase in observed disc ra- 
dius if this is tracing where the disc has recently been stirred. While 
in principle this would provide a way to discriminate between dif- 
ferent stirring mechanisms (for self-stirring r oc t 1 ' 3 while for 
planet-stirring r oc t 2 ' 9 ), the difference is so small that it would 
be very difficult to distinguish in practice. In any event, there is 
not currently evidence fo r any dependence of disc radius on time 
dNaiita & Williamsl2005h . 

6.2.2 Discs unlikely to be self-stirred 

Some stars are young (~ 10 Myr) yet already have bright de- 
bris discs of large radius (~ 100 AU). Such discs cannot have 
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self-stirred unless the disc is sufficiently dense, and in such sys- 
tems planet-stirring may be a viable alternative. To quantify this, 
we calculated the minimum x m required for a disc to self stir in 
less than the system age, for discs around 23 FGK sta rs and 35 A 
stars with published 24/25 and 70/60 micron excesses l Moor et all 
20061: IBeichman et alj l2006t ISu et al.1 l2006t iTrilling et al.l 120071 : 
Hillenbrand et al J2008h . Disc radii were estimated by fitting black- 
body curves to the IR excess. For FGK stars we increased these 
radii by a factor of three because a comparison with the radii known 
directly from those discs which have been imaged showed that the 
black body fits systematically underestimate the radii by roughly 
this amount; this is likely due to the small blow-out size for dust 
in these discs. If radii were available from imaging, we used these 
in preference to the black-body fits. We identified that a disc would 
have trouble self-stirring if a; m , m in > 10. 

Very massive discs would also have been gravitationally un- 
stable when gas was still present. We can calculate the minimum 
density for gravitational instability using the Toomre criterion Q = 
n C g™ < 1 for instability, where n is the mean motion of the disc, 
c s is the sound speed in the disc gas, and S E is the surfa ce density 
of gas. The sound speed c s ~ (h/a)v-^ cp (Pringle 1981), where h 
is the disc scale height. Assuming h/a = 0.1 and a dust:gas ratio 
of 1:100, so that S g = 100E with £ given by Equation l|29]l, this 
gives 

-1/2 



anqular distance (arcsec) 
0.1 1.0 



Q ^ 10' 



(l Au) 



(31) 



for instability. For a disc at 100 AU, this corresponds to a maximum 
i m of ~ 10. Higher Xm would still be possible through metallicity 
enhancement without affecting the gas mass and therefore gravita- 
tional stability. 

We identify two discs with a minimum surface density for self- 
stirring a; m , min ^ 10: HD 181327 (a; m ,min = 17) and HD 202917 
(£m,mm = 10). Both these disc s have been ima ged, with r adii 
a = 86 AU dSchneider et al.120061) and a sa 80 AU Jkristl2007l) re- 
spectively. It may be that such discs do indeed have a; m , m i n 10; 
i.e., they may be at the top end of the disc mass distribution, in 
which case they may be self-stirred, assuming that they have man- 
aged to avoid the gravitational instability mentioned in the previous 
paragraph. In the case of HD 181327, however, there is independent 
evidence in support of plane t stirring: the disc has an azimuthal 
asymmetry l ichen et alj2008l) which could be due to planetary sec- 
ular perturbations dWvatt et alj|l999h . so this system in particular 
warrants further investigation. In such a disc, we can place con- 
straints on the parameters of the planet responsible for stirring by 



requiring t c 



< fa 



. Figure [8] shows that such planets are likely 



to be be far from the star. This, together with the host stars' youth, 
makes them good targets for direct imaging. 



6.3 Exoplanet population 

We now attempt to ascertain whether there is any observational evi- 
dence for planet stirring. We begin by looking at a statistical sample 
of exoplanets. If planets are a common cause of disc stirring then 
we might expect there to be a correlation between the planetary 
parameters $ and a* and the presence of infrared excess, higher 
values of these parameters correlating with IR excess. Of the two, 
a* is the more fundamental because it describes the planet's ab- 
solute ability to stir a disc within the context of the planet-stirring 
model, independently of any other sources of stirring. A star host- 
ing a planet may also host a planetesimal disc. If it does, and if 
planet-stirring were the sole stirring mechanism, then we would 
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Figure 8. The shaded region shows the region of a p i-m, p i parameter space 
a perturbing planet must occupy to stir the disc in the HD 181327 system 
in less than the system age of ~12Myr, assuming a planetary eccentricity 
of 0.1 or 0.5. Planets with higher eccentricity can stir the disc if they have 
lower mass or semi-major axis than ones with lower eccentricity. The disc 
is indicated with the thick vertical line. 



expect only those planets with high enough a* to exhibit IR excess. 
The parameter $ quantifies the relative importance of the planet- 
and self-stirring models, so the interpretation of any correlation be- 
tween high $ and a disc, should one exist, is not so clear. 

We take 57 planet-hostin g stars w ith published Spitze r /MIPS 
photometry (IBeichman et al. I [20061: iMoro-Martm et"ai] |2007|; 
ITrilling et al.ll2008l ; iBackman et al.l 120091) . Our sample is identi- 
cal to that of Bryden et al. (in prep.), but with the two M dwarfs 
GJ 436 and GJ 876, and the G dwarf HD 3363 6 whose compan - 
ion has been determined to be of stellar mass jBean et al.li2007i> . 
removed. Of these 57 stars, 10 show significant excess emission at 
24 fim a nd/or 70 /xm, and ar e classed as disc hosts. Exoplanet data 
are from Butler et al] d2006t) for most planets, except fo r e Eridani 
jBenedict et ai]|2006l) and HD 69830 dLovis et al.ll2006l) . 

In Figure l[9} we plot parameters $ and a* for the 57 systems. 
For multi-planet systems, we plot <E> and a* for the planet which 
stirs the disc quickest, treating the system as if the planets' pertur- 
bations acted independently; see i]6.4.3l for a more detailed discus- 
sion. We also plot the line for constant e p i = 0.1, m p i = 10 -3 Mq 
(a 'typical planet'), and varying a p i, which fits the points rather 
well, showing that most of the range of $ and a* can be accounted 
for by spread in a p \. The planet's semi-major axis is the most im- 
portant parameter in determining whether a disc can be self- or 
planet-stirred. The region on the right shows the region of param- 
eter space in which a planet can stir an a = 10 AU, x m — 1 disc, 
according to Equations J3 0| > and l |26t . 

Figure l[9} does not suggest a difference in the distributions of 
$ and a* between planets orbiting disc hosting and non disc host- 
ing stars. This is confirmed by Kolmogorov-Smirnov tests. The p- 
value j3 from one-dimensional KS tests comparing the disc hosting 



J The p- value gives the probability of observing a more extreme test statis- 
tic under the assumption of the null hypothesis: that the distribution of the 
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Figure 9. Plot of a* against <E> for systems in the sample of Bryden et al. (in 
prep.)- We plot only the planet which stirs the disc quickest for multi-planet 
systems. Black crosses show planets orbiting stars without discs. Blue stars 
show planets orbiting stars with discs. Dotted line marks 10 AU. Dashed 



line marks a* against <E> for planets with e p j 
and varying a p \. We set ic m = 1. 
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and non disc hosting samples are 0.986 when comparing the dis- 
tributions of <J> and 0.917 when comparing the distributions of a* . 
Thus the distributions of planetary parameters are statistically in- 
distinguishable between disc hosts and non disc hosts. 

The implications of this for the relative importance of self- 
and planet-stirring are however unclear. Due to the many processes 
doubtless taking place, it is likely that should any correlation be 
present it has been diluted. Larger samples, at a range of ages (all 
but two of the stars in this sample are over 1 Gyr old), may be nec- 
essary to properly determine the evolutionary processes at work. 

6.4 Case studies 

We now proceed to examine some individual systems, categoris- 
ing them somewhat arbitrarily by the planet's semi-major axis. We 
conclude this subsection by briefly looking at multi-planet systems. 

6.4.1 Jupiter analogues 

e Eridani hosts a highly eccentric (e P i = 0.7) Jupiter-m ass (m p i = 

1.5 x 1O" 3 M ) planet at 3.4 AU dBenedict et alj2006f) . This planet 
can stir any disc of 80 m planetesimals out to a* = 3000 AU. 
e Eridani also hosts a cold debris disc extending from 35-1 10 AU, 
with surface brightness peaking at 60 AU, as wel l as at least one 
unresolved warm inner belt dGreaveset all 120051: 1 Backm an et al.l 
2009). Clearly, all the dust is well within the maximum radius for 
planet-stirring by the a* criterion. We can also compare t C ross to the 
star's age, estimated at ~ 850 Myr dDi Folco et al.ll2004[) . Because 
Across ~ 40 Myr at 60 AU, we see that planet-stirring by planet b 
is inevitable within the system age. Figure dlOt shows the region 



parameter is the same for both samples. Small p-values suggest that the 
populations have different distributions. 



Figure 10. Planet parameters needed to stir e Eridani's disc at 60 AU, within 
the age of the system. The diagonal lines show the minimum planet mass 
at each semi-major axis required to stir the disc for e p i = 0.1 and e p i = 
0.7. We also show detection limits for RV (dotted line) and direct imaging 
(dashed line). The real location of e Eridani b is marked with a cross. 



of parameter space within which a planet must lie if it is to have 
stirred the disc within 850 Myr, for both a fiducial planetary eccen- 
tricity e p i = 0.1 and the real planetary eccentricity e p i = 0.7. The 
figure also shows the parameter space accessible to 20 years' radial 
velocity observations at 15 m s _1 precision, and the rough sensitiv- 
ity l imits from the direc t imaging searches of lMarengo et al.l ( 2006) 
and |jansonetalJ l l2008l) . 

We can also calculate the maximum size of bodies that can be 
destroyed by planet stirring by e Eridani b. At 60 AU, Equation i25l 
gives J?max ~ 110 km. The planetary perturbations therefore have 
an impact over a wide range of the size distribution of planetesi- 
mals. 

It is intriguing that the planet e Eridani b is only just able to 
stir the disc at 110 AU, within 850 Myr. This may be coincidence, 
but may hint that there is an as yet unstirred disc region beyond 
110AU, which the planet's perturbations have not yet reached. It 
is, however, worth noting that the disc could also be stirred by any 
other planet in the syst em, such as that po stulated to explain the 
disc's clumpy structure dGreaves et al.l2 005), although here the in- 
teractions would not be purely secular. 

6.4.2 Neptune analogues 

Fomalhaut' s lon g-suspected planet was recently imaged by 
iKalas et al] d2008l) . Its orbital elements are estimated at e p \ = 
0.11, and a D i = 115 AU , with an upper limit of 3Mj for the 
mass dChiang et al. |l2009l). The disc lies at a radius of ~ 140 AU 
dKalas, Graham & Clampinl2005l) . At this radius, the time-scale for 
planet stirring is only 0.65 Myr, orders of magnitude les s than the 
star's age (~ 200 Myr, iBarrado v Navascues et alj|l997l) . We also 
find a* = 1.2 x 10 4 AU, so the disc is well within the radial limits 
for planet-stirring. 

Figure Jilt shows the planet masses and semi-major axes re- 
quired to stir the Fomalhaut disc within the system's age, assuming 
a planetary eccentricity of 0.1. We can see that even if the planet's 
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Figure 11. Planet parameters needed to stir the disc of Fomalhaut. The 
shaded region marks the planet parameters where disc stirring can occur 
within the system age. The dashed line marks the planet parameters which 
can stir the disc as quickly as Fom b (0.65 Myr). The planet's semi-major 
axis and maximum mass are marked with an arrow. 



mass is significantly less than the maximum of ~ 3Mj, planet- 
stirring would occur within the age of the system. 

However, Fomalhaut presents two complications. Firstly, there 
are difficulties with in situ formation of the planet Fom b because 
the time-scales for core accretion are so long: recall that the self- 
stirring time-scale gives the time required to form Pluto-sized ob- 
jects in situ, while the mass of Fom b may be as high as that of 
Jupiter. The time to form a Pluto-sized body at Fomalhaut b's or- 
bit is around 150 Myr for x m — 1. The planet most likely formed 
closer to the star and later mo ved to its current lo cation, for example 
by outwards migration (e.g., Martin et alj|2007t) or being scattered 
by another planet (e.g., Veras. Crepp & For j|2009l) . Both of these 
processes would however likely disturb the disc as well. 

Secondly, we note that, although et ~ 0.1, the material in 
the Fomalhaut disc appears to hav e very low proper eccentricities 
( Guillen 2006; IChiang et alj |2009), as evidenced by the sharp in- 
ner edge to the disc. If the proper eccentricity of Fomalhaut's disc 
is only 10 per cent of the forced eccentricity then this increases 
the time-scale for orbit crossing to 6.5 Myr (see Equationll4t. still 
much less than the system's age. Reducing proper eccentricities 
also reduces the relative velocities amongst planetesimals in direct 
proportion, although given the large value of a* this will not pre- 
vent Fom b from causing erosive collisions. 

For the Solar System's Neptune we find a* = 730 AU, mak- 
ing the Kuiper Belt able to be stirred by Neptune. However, when 
we compare with self-stirring we find that planet-stirring acts more 
quickly only out to $ = 33 AU, so Neptune's secular perturbations 
would not have stirred the belt before Pluto formed, assuming that 
the planets formed at their current semi-major axefl So this sim- 
ple model is consistent with the outer Solar System, although we 
note that the dynamical evolution of the early Kuiper Belt and outer 



4 We find similar values for other giant planets; e.g., for Jupiter, a* 
720 AU and* = 22 AU 



Figure 12. Solid lines: Precession rate A for planetesimals in the Sun- 
Jupiter-Saturn system. The horizontal lines mark the eigenfrequencies. Dot- 
ted lines: precession rate and eigenfrequencies for the same system but with 
Saturn's mass reduced to that of Earth. Dashed line: precession rate for 
planetesimals perturbed by Jupiter alone. Secular resonances occur when 
the precession rate equals one of the system's eigenfrequencies; there are 
no secular resonances in the case of a single perturber. 



planets may have been more complicated than in-sit u formation of 
Neptu ne followed by growth of Kuiper Belt Objects iTsiga nTs et al.l 
2005). We also note that highly excited eccentricities and inclina- 
tions of KBOs may have been re quired to explain the details of th e 
capture of Neptune's Trojans dNesvornv & VokrouhlickH l2009h . 
and capture of KBO s into high order mean motion resonances (e.g., 
IChiang et alj|2003l) . Such high inclinations might be achievable 
through self-stirring but not planet-stirring. 



6.4.3 Multi-planet systems 

When dealing with multiple planets previously we treated the disc 
as being stirred by the planet with the lowest t cr0 aa, assuming that 
the other planets had no effect on the disc. Such an approach is 
unrealistic because it neglects not only the effects of other planets 
on the disc, but also the mutual interactions of the planets amongst 
themselves. 

We plot the precession rate A for planetesimals orbiting in the 
Sun- Jupiter-Saturn system in Figure d!2t . This also shows the lo- 
cation of secular resonances, where the planetesimal's precession 
rate equals one of the system's eigenfrequencies and the forced ec- 
centricity is formally infinite. 

Figure U21 also shows the effect of reducing Saturn's mass 
to that of Earth: the precession rate approaches that in the single- 
planet case of Jupiter alone, and the width of the region strongly af- 
fected by the outer planet's perturbations decreases. So as far as the 
precession rate is concerned, the behaviour is similar to the single- 
planet case. 

Performing a similar analysis to that in Sj3] we find that, for 
planetesimals on initially circular orbits, the time-scale for orbit 
crossing in the multi-planet case is given by 
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Figure 13. Solid line: orbit-crossing time-scale (cross for planetesimals in 
the Sun-Jupiter-Saturn system. Dotted line: same, with Saturn's mass re- 
duced to that of Earth. Dashed lines: t cr oss assuming a one-planet system 
with Jupiter stirring the disc. Dot-dash line: the asymptotic approximation 
to the solid for small a. As well as being short near the planets (marked with 
their initials), t CTOBB is also small near the four secular resonances (marked 
with vertical lines for the Sun-Jupiter-Saturn system). 
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where bj — 7/2 or —3/2 when planet j is an internal or external 
perturber respectively. 

The stirring time t cr0 ss for the Sun- Jupiter-Saturn system is 
plotted in Figure l |13t . We also show t C ross with Saturn's mass re- 
duced to that of Earth, and for the single-planet case with Jupiter 
alone perturbing the disc. With Saturn at its true mass, we see that 
the crossing timescale is greatly reduced close to the planets. How- 
ever, beyond the outermost secular resonance, the dependence of 
t cr0 ss on a steepens. Specifically, for large a, t CIOBB oc a 8 rather 
than a 4 ' 5 . This is because we now have \A\ < \gi\ in the forced 
eccentricity term in Equation i32\ . and so the a dependence of A 
and Bj no longer cancels. For planetesimals beyond 20 AU, in- 
troducing another perturber has increased the time-scale for orbit 
crossing. When Saturn's mass is reduced to that of Earth, we see a 
large region between Saturn and the outer secular resonance where 
the the time-scale is the same as for the case with Jupiter alone: 
because of the large disparity in masses, the perturbations are dom- 
inated by Jupiter. 

As an example of a multi-planet system with a debris disc, 
consider HD 38529. This star hosts a 0.8 Jupiter mass planet on 
a 0.13 AU orbit and a 12.2 Jupiter mass planet on a 3.74 AU or- 
bit. The secular dynamics of the system, incl uding both planets 
and m assless planetesimals, were modelled by iMoro-Martm et alj 
(2007), who concluded from the dynamical analysis and SED fit- 
ting that the planetesimals reside in a dynamically stable region 
at 20-50 AU between secular resonances. Figure d 1 4b shows the 
crossing time-scale for HD 38529. Within the region 20-50 AU, 
the crossing time-scale is close to that achieved by planet c alone, 
due to its higher mass and larger semi-major axis. Within this re- 
gion of the disc, the planets' secular perturbations induce crossing 



Figure 14. Solid line: orbit-crossing time-scale tcross for the HD 38529 
system. Planets are shown with their letters, and the secular resonances are 
marked with vertical lines. The dashed line shows the crossing time if only 
planet c is stirring the disc. The sharp peak at ss 0.7AU occurs where 
dA/da = (see text). 



of neighbouring initially circular orbits on time-scales of < 1 Myr. 
It may well be the case then that there are no bodies larger than 
several kilometers in radius in the disc (but see i]6.5.2t . 

Note that t C ross is formally infinite when dA/da = 0. This 
would appear to mean that there is a particular semi-major axis be- 
tween the planets where the perturbations can never induce orbit- 
crossing. However, this singularity is merely a mathematical arte- 
fact (see ij3}: in reality, this region can still be stirred by secular 
perturbations, given sufficient time. 



6.5 Limitations of model 

We have presented a simple picture of the effect of secular pertur- 
bations on a planetesimal disc. Here we clarify the assumptions and 
limitations attached to our model. 



6.5.1 Collision model 

To determine the outcome of collisions, we have chosen one partic- 
ular scaling law for threshold collision energy. Other scaling laws 
differ in the planetesimal radius at which the minimum of is 
attained, and the value of the minimum itself. Given the very large 
values of a* for 100 m bodies, this is unlikely to be important for 
this size of planetesimal, but may be important if the planetesimals 
are larger. 

We have assumed that collisions occur as soon as the orbits be- 
gin to cross, and accounting for collision rates will slightly increase 
the stirring time. 



6.5.2 Initial conditions 



In common with other studies (e.g., Theb ault et al] 120061 : 
iKenvon & Bromlevl [2008), we have taken the initial planetesi- 
mal orbits to be circular. Despite promising recent progress (e.g., 
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Johans en et al ]|200l l Cuzzi et al]|2008h . the formation of planetes- 
imals is still not fully understood, and so these initial conditions, 
although reasonable, are not rigorously justified. In particular, we 
note that planetesimals can acquire moderate (~ 0.05) eccentrici- 
ties if they orbit within a protoplan etary disc that has density fluctu- 
ations induced b y its self-grav ity teritsch. Clarke & Lodatoll2008t) 
or by the MRI (Nelson 200^). Starting at different eccentricities 
can be dealt with by changing the proper eccentricity e p in Equa- 
tion d!4t . so long as e p is a single-valued function of semi-major 
axis. 

We have also introduced the planet instantaneously, assuming 
that it forms at time t = 0. Statements about the stirring time-scale 
should therefore be qualified by including the time taken for the 
planet to form, during which time the planetesimals will themselves 
be growing. This will affect the maximum radius of bodies which 
exist in the disc when the orbits begin to cross. 

6.5.3 Unmodelled processes 

Dynamically, we have neglected any non-secular dynamical effects 
of the planet on the disc. While this is valid for the planetesimals far 
from the planet, when they are close they begin to experience reso- 
nant interactions as well as secular (see Figure^. Furthermore, the 
asymptotic expression for crossing time-scale (Equation 1 151) is no 
longer valid. In fact, this approximation over-estimates the time- 
scale because, as a — > a P i, we find A — » oo while ei — > e p i 
dMurrav & Dermott| [T999). and so t C ross — > 0, whereas extrapola- 
tion of the asymptotic expression predicts a finite value. So both 
simplifying the secular interactions and neglecting non-secular in- 
teractions will tend to overestimate the crossing time-scale. Thus 
the simplified dynamics underestimates the ability of planets to stir 
discs located close to them. This may be particularly relevant for 
such sys tems as Fomalhaut, where several strong resonances lie in 
the disc JChiang et alj2009h . 

We have neglected any damping in the disc. Collisional damp- 
ing has been inv oked to explain the low proper eccentricit ies in 
Fomalhaut's disc jOuillen & Fabeill2006l : Ichiang et alj|2009h . Re- 
ducing proper eccentricity increases both the time-scale for orbit 
crossing and the relative velocities, as we have already described. 

Finally, while we have focussed on planetary secular pertur- 
bations as a disc stirring mechanism, and compared them to in situ 
planet formation, it is important to realise that there may b e other 
causes of dynamical excitation. iKenvon & Bromlevl {2002) inves- 
tigated the effects of a stellar flyby, but found that the perturba- 
tions were rapidly damped by collisions. Other mechanisms which 
have not been so thoroughly investigated in the context of debris 
disc evolution include planet formation proce eding more rapidly 
when a gap-ope ning giant planet has formed IBrvden et al . 2000; 
iThommej 12005). or an outer disc being stirred by planetesimals 
that ha ve been scattered out from the inner system ( Goldreich et al. 
l2004ah . 



7 SUMMARY AND CONCLUSIONS 

Our main conclusion is that a planetesimal belt at several 10's of 
AU can be stirred by an eccentric giant planet at only a few AU. 
Debris discs do not require any bodies larger than a few km in size 
beyond a few AU from the star to explain the observed dust pro- 
duction. 

To reach this conclusion, we investigated the effects of sec- 
ular perturbations from an eccentric planet on a dynamically cold 



disc, to assess whether this might be the origin of debris disc stir- 
ring. Over time, neighbouring orbits acquire sufficiently different 
eccentricities and longitudes of periapse that they begin to cross. 
We derived an analytical expression for the time t cr0 ss for this to 
occur (Equation 1 16b. which agrees well with previously published 
N-body simulations. After this time, the planet's secular perturba- 
tions quickly impose a mean relative velocity (v re i) ~ 1.4efi>ke P . 
This is similar to that for a planetesimal swarm without external 
perturbers, despite the high degree of apsidal alignment forced on 
the orbits by secular perturbations. 

When the relative velocities increase, the disc may brighten as 
a result of increased dust production, and further growth of plan- 
etesimals may be inhibited or halted. We derived an expression for 
a*, the maximum range of semi-major axes over which a planet's 
perturbations can destroy planetesimals (Equation 1 26b. The range 
a* increases with planetary eccentricity and strongly increases with 
planetary semi-major axis. It is also a function of planetesimal size. 
For the weakest planetesimals (~ 100 m in radius), we found typi- 
cal values of a* of several hundred AU. 

We then compared the time-scales for pl anet stirring with the 
time-s cales for self stirring from the models of lKenvon & Bromlevl 
(2008). Because t C ross has a stronger dependence on a than does 
tp\, assuming that only these two proceeses are operating, we found 
that typically the disc closer to the planet will be planet-stirred, and 
the disc further away will be self-stirred, and we identified another 
parameter $ which demarcates the outer reach of the planet's per- 
turbations, beyond which the self-stirring time-scale is shorter. This 
parameter $ is typically much smaller than a*, and whether a disc 
can be planet-stirred before it is self-stirred depends on the disc 
density as well as planet parameters. 

For a sample of RV planet hosts observed by Spitzer, we find 
no correlation between the magnitude of $ or a* and the presence 
of a disc. While this may suggest that planet stirring is not ubiqui- 
tous, the degree of correlation we would expect is not clear. 

However, for some individual systems it appears that a known 
planet will have stirred the disc on a time-scale shorter than the sys- 
tem age and/or before it it likely to self-stir. In particular, we iden- 
tify e Eridani, Fomalhaut, and HD 38529 as being in this category. 
Future studies of such discs should take into account the effects of 
planetary perturbations on the disc's evolution, collisional as well 
as dynamical. We also speculate that the discs of HD 181327 and 
HD 202917 may have been stirred by as yet undiscovered planets, 
since the only other stirring mechanism proposed would only work 
for discs with high surface densities 10 MMSN). 

Future work in this area should further investigate the effects 
of planetary perturbations on the collisional evolution of a disc, 
with a view to (a) clarifying the extent to which the perturbations 
can inhibit the further growth of planetesimals and (b) determining 
the evolution of dust production and hence IR luminosity. The lat- 
ter in particular will enable valuable observational tests of planet- 
stirring and its role vis-a-vis other stirring mechanisms. We also 
note that, while we have focused on internal perturbers in our dis- 
cussion, our formulae for external perturbers will be relevant for 
investigations into the effects of binary companions on discs. 
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